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Abstract. Teleportation of optical field states (as continuous 
quantum variables) is usually described in terms of Wigner func- 
tions. This is in marked contrast to the theoretical treatment 
of teleportation of qubits. In this paper we show that by using 
the holomorphic representation of the canonical commutation rela- 
tions, teleportation of continuous quantum variables can be treated 
in complete analogy to the case of teleportation of qubits. In order 
to emphasize this analogy, short descriptions of the basic experi- 
mental schemes both for teleportation of qubits and of continuous 
variables are included. We conclude our paper with a brief discus- 
sion of the effectiveness of our description of continuous variable 
teleportation and of the role of localization of quantum states in 
teleportation problems. 



1. Introduction 

The essence of the quantum information processing consists in (1) 
generating quantum entanglements among quantum systems, (2) con- 
trolling the quantum entanglements. Quantum teleportation contains 
the essence of quantum information processing ((1) and (2)). Quan- 
tum teleportation will play an important role in the realization of the 
quantum computer. 

Quantum teleportation is first discovered in 1993 by Bennett et al. 
[2]. A few years later two experimental reports on quantum teleporta- 
tion were published, in 1997 by Bouwmeester et al. [I], and in 1998 by 
Furusawa et al. [ID] . 

The experiment (by Bouwmeester et al. [I]) for the teleportation of 
photon states (qubits) is difficult, because the efficiency of single photon 
experiments is presently restricted in principle due to the inability to 
identify all four Bell states, and also in practice by the low efficiency 
of single photon production and detection. In contrast, the important 
feature of the technique used in the experiment for the teleportation 
of optical field states by Furusawa et al. [10] is its high efficiency. 
This is due to the in principle ability to perform the required joint 
measurements and the technical maturity of optical field detection. 
The experiment of Furusawa et al. [10] is often considered to be the 
first experimental realization of quantum teleportation. 



2 



S. NAGAMACHI AND E. BRUNING 



The purpose of this lecture is to present teleportation of photon 
states and optical field states in a unified way, with emphasis on a new 
approach to teleportation of optical field states (quantum teleportation 
of continuous quantum variables), based on the holomorhic represen- 
tation of the canonical commutation relations instead on the use of 
Wigner functions. Here we try to have our lecture self contained as far 
as possible (not with respect to the underlying literature but respect 
to the arguments which are used). 

The mathematics for the teleportation of photon states (qubits) is 
the theory of 8 x 8 matrices, that is, the theory of 8 dimensional Hilbert 
spaces. However, in order to describe the optical field states we need the 
theory of infinite dimensional Hilbert spaces. The experiment uses the 
entangled states of squeezed laser beams with squeezing parameter r. 
Since we can neither generate an infinitely squeezed (r = oo) EPR state 
nor prepare ideal detectors with efficiency 1, we cannot have complete 
teleportation tp out = ip- m . We have to measure the quality of the output 
state Vout- To do so, one often uses the notion of fidelity which is 
defined by F(p, a) = tr a/ p-^Vp 1 / 2 for density operators p and a. 

In the description of teleportation of continuous quantum variables, 
one often uses the Wigner function which allows a compact form for 
the output state ip out for a given coherent input state ip- m and has 
an intimate connection with the fidelity (see [U]). Such a formulation 
sets the theoretical description teleportation of continuous quantum 
variables apart from the description of teleportation of qubits which is 
based on the use of state vectors and projective measurements. 

In order to understand the essence of the teleportation of continu- 
ous quantum variables and to show the close analogy to teleportation 
of quibits, we use the "holomorphic representation" of the canonical 
commutation relations (CCR), which was introduced by Bargmann 
[lj. This representation allows to derive easily explicit formulae for 
the theoretical description of the basic operations and objects used in 
the experiment. 

The contents of this lecture can briefly be described as follows: In 
Section 2, we recall quantum teleportation of qubits, in Section 3, we 
present briefly the experiment of teleportation of qubits. In Section 
4, we introduce the holomorphic representation of CCR and deter- 
mine explicitly the kernels of the mathematical operations (Bogoliubov 
transformation) which are needed in the Section 5 to derive the mathe- 
matical realization of the devices for the manipulation of photon states 
(squeezed vacuum state, half-beam splitter, displacement operator). 

The theory of laser and parametric oscillator (amplifier) are impor- 
tant for quantum teleportation but not included in this lecture. We 
refer for example to p3J IZj- In Section 6, we present our approach 
to quantum teleportation of continuous variables. In Section 7, the 
experiment made by Furusawa et al. [10J is briefly presented with a 
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discussion of a controversy about this experiment. Some remarks on 
the notion of locality in the theory of quantum teleportation conclude 
this lecture. 

For the convenience of the reader, the first part of an appendix ex- 
plains the notion of generalized states (which play an essential role in 
our approach) and mentions some basic properties; the second part 
gives the detailed proofs for the results presented in Section 4. 

2. Teleportation of qubits 
For the teleportation of qubits 

M = a\0) + (3\1) # ( ° ) G C 2 , |0> = ( I ) , |1> = ( J ) 

the following single qubit gates are used: 

x =(° ;).*(;)=(f).*(io).ii>)=(n>.io». 

""(? 1)'*-( 1 -1).'-7S(I -1 ) = <*+*>/ A 

ff(|0),|l» = -^(|0> + |l>,|0>-|l». 

The matrix H is called the Hadamard gate. 

In order to treat multi-qubits, we must consider composite systems. 
Axiom of composite system: The state space of the composite phys- 
ical system is the tensor product of the state spaces of the component 
physical systems. Moreover, if we have systems numbered 1 through 
n, and system number i is prepared in the state \ipi), then the state of 
the total system is <g> \ip2) <S> • • • <S> \ipn)- 

The two qubit gate (Controlled-Not gate) Mcnot will play an impor- 
tant role in our lecture. This gate acts on the basis vectors 

\ij) = \i) ® \ j) = 

as follows: 

M CNOT (|00), |01), 1 10), |11)) = (|00), |01), |11), |10». 

2.1. EPR pair, Bell states. The essence of quantum information 
science is quantum entanglement and its manipulation. The entan- 
gled states called EPR pair or Bell states are created by using 
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Hadamard gate 77 and controlled- not gate Mcnot : 
(2.1) 

M CNOT (77 ® 7)|00) = -^=M C not(|0) + |1» ® |0) = ^=(|00) + |11» = |/3 o), 
M CNOT (H® J)j01> = -Lm cnot (|0) + |1» ® |1) = -^=(|01> + |10>) = \(3 0l ), 
M CNOT (#® /)|10) = -^M CNOT (|0) - |1)) ® |0> = ^=(|00) - |11)) = |Ao), 

M CNOT (# ® 7)|11) = -Lm cnot (|0) - |1)) ® |1) = -i=(|01> - |10» = |/3 n >. 

Remark 2.1. The canonical basis can be written in terms of the 
Bell basis |/%) as follows: 

|00> = -^=(|/?oo> + |/3io», |H) = ^=(|/?oo) - |Ao)), 
(2.2) f f 

|01) = -j=(\(3 01 ) + |/3 n », |10) = -/=(|A>i) - 

2.2. Description of quantum measurement. Axiom of quan- 
tum measurement: Quantum measurements are described by a col- 
lection {M m } of projection operators which appear in the spectral de- 
composition of the observable 



M = ^ mM n 



The index m refers to the measurement outcomes that may occur in 
the experiment. If the state of the quantum system is \ip) (HIV') II = 1) 
immediately before the measurement then the probability that result 
m occurs is given by 

p(m) = (ij\M m \ij), 
and the state of the system after such an ideal measurement is 



The projection operators satisfy the completeness equation, 
M m = I = Identity operator on the state space. 

m 

2.3. Quantum teleportation. Now we can describe the process of 
quantum teleportation, which is illustrated in Fig. 1. 

(1) Alice and Bob prepare an EPR pair \/3oo)ab 
|A)o) ab = ^(\0)a ® |0)b + |1)a ® \1)b) = ^=(\0)a\0)b + \1)a\1)b) 

by applying a Hadamard gate and a controlled-not gate to the 
state |0)^4 <g> |0) B , the composite state of Alice's state \0)a and 
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Figure 1. Quantum circuit of teleportation. Victor's 
state \ip)v is given to Alice and reproduced at Bob's lab- 
oratory. 



Bob's state |0) B . Now Alice and Bob share the state \Pqo)ab 
(|0)a, 1 1)^4 are Alice's states, and \0)b, are Bob's states). 
Victor gives Alice a state \ip)v (IHV')vll — 1) to send to Bob. 
The state of the total system is 

l^o) = W)v ® \Poo)ab = (a\0) v + p\l) v ) ® -±=(\0) A \0) B + \1)a\1)b). 

The state \ipo) = \^)v ® \Poq)ab can be rewitten as: 

l^o) = (a\0) v + P\l) v ) ® ^=(\0) A \0) B + \1)a\1)b) 

= ^=[a\0) v \0)A\0)B + a\0) v \l) a\1)b 
+/3|1)v|0) a |0) b + /3|1)v|1)a|1)b] 
= \h\Poo)va + \Pio)va)\0)b + a(\Poi)vA + \Pu)va)\1)b 
+P(\Poi)va-\Pii)va)\0)b + P(\Poo)va-\Pio)va)\1)b] 
= ^[\(3 00 )vA(a\0) B + P\1) B ) + \Poi)va(<x\1) b + P\0)b) 

(2.3) + \(3 10 } VA (a\0) B - P\1) B ) + \Pn)vA(a\l) B - P\0) B )\. 
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(2) Alice performs the Bell-state measurement, a measurement which 
determines |/3y), and the state of the system after the measure- 
ment is 

\(3oo)vA(a\0) B + (3\l) B ) if (ij) = (00) 
|Au)™(a|l>B + /3|0)fl) if(u') = (01) 
\P W )vA(a\0)B-P\l) B ) if(u) = (10) ' 
\l3n)vA(a\l} B -l3\0}B) if (ij) = (H) 

(3) Alice sends the classical information (ij) to Bob. Then Bob 
sends Bob's qubit through I, X, Z, XZ according to the result 
(00), (01), (10), (11), obtaining a\0) B + (3\1) B . 

(4) Since the operator (if® 7)Mcnot is the inverse of the operator 
Mcnot (H®I), (if <g) I)M CNOT sends the Bell basis \Pij)vA to 
the canonical basis \ij)vA = \i)v® \])a of product states. Thus 
Alice's Bell-state measurement is performed by sending the Bell 
basis to the canonical basis \ij)vA and making the measurement 
determining \ij)vA- In practice, Alice sends her qubits through 
a CNOT gate, and then the first qubit through a Hadamard 
gate if, obtaining 

^[\Q0)vaH0) b + (3\1)b) + \01) VA (a\l) B + (3\0)b) 

(2.4) + |10)™(a|0) fl - f3\l) B ) + |ll) VA (a|l) fl - /3\0) B )]. 

Let Mij = \ij)vAVA{ij\ ® I B , (hj — 1,2). Then these projec- 
tion operators satisfy the completeness equation. 

(5) Alice performs a measurement {M^}, i.e., measures the observ- 
able M = E! j= o ( 2i + j) M ij- The probability that result (ij) 
occurs (equivalently m = 2i + j) is 

P (ij) = (HMijIfa) = \, 

and the state of the system after the measurement is 

f |00) VA (a|0)B + /3|l>B) if (ij) = (00) m = 

MijM I |01)vA(a|l>B + )9|0)B) if(u) = (01) m=l 

JWMM I |10> V A(a|0>B-/3|l>fl) if (ij) = (10) m = 2 ' 

[ |ll) VA (a|l)fl-)9|0>fl) if (zj) = (11) m = 3 

3. Experiment of the qubit teleportation 

In this section we discuss the experiment performed by Bouwmeester 
et al. in 1997 [I] using qubit states. Two pairs of entangled photons 
are generated by a polarized non-degenerate parametric process. Let 
|0) represent the horizontally polarized single photon state | <->) and 
|1) the vertically polarized single photon state | £). 
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For these polarization states we calculate the Bell basis according to 
( 12. ip . Then the state \ipo) = \ip)v^> \Poo)ab is rewritten as follows, (see 

my- 

l^o) = 1[\Mva(u\0) b + p\l) B ) + \I3oi)va(oc\1)b + P\0)b) 

+ \Pio)vA(a\0) B - 0\1)b + \Pn)vA(a\l) B - P\0) B )]- 

A Bell-state measurement is possible for the state \/3u)va, under the 

condition that there are at most one photon on each mode of Victor and 

Alice. Let ayj, ctAk be annihilation operators of polarization j, k =<->, \. 

The half-beam splitter causes the Bogoliubov transformation (see f)5.4p ) 

1 1 
b j = -j=(a V j + a Aj ), by = -j=(—a V j + aAj) 

Let |0) be the vacuum for ayj, a Ak . Then we have 



1 



(Q\b lk b 0j = -(Q\(-a vk + a Ak )(a Vj + a A 



j ■ 



We assume that there is at most one photon in each mode (Victor's 
mode or Alice's mode). On such a condition, we can perform the 
Bell-state measurement by the simultaneous photon counting after the 
half-beam splitter. For such a state we can ignore ayk^vj and 
a>Ak&Aj we have 

(Q\b lk b 0j \^) = (Vt\(-a Aj a Vk + a Ak a Vj )\ip) = ±-^= VA (f3 11 \ip) 
for j ^ k, and (fl\bi k b j\ip) = for j = k, where we used the fact 

(- a Vk a Aj + a Vj a Ak)\ n ) = ±^\Pll)vA- 

This shows that the simultaneous photon detection is equivalent to 
va^hI) that is, equivalent to the Bell-state measurement. However, 
this measurement cannot identify the other three Bell states, and has 
a fatal drawback. The photon counting technique of today cannot dis- 
tinguish whether only one photon is coming or more than two photons 
are coming simultaneously. So, if two photons come in Victor's mode 
simultaneously, we cannot neglect the term ay k avj, the measurement 
is not the Bell-state measurement. 

4. Holomorhic Representation of CCR 

The "holomorphic representation" of the canonical commutation re- 
lations (CCR) was introduced by Bargmann [lj for the finite dimen- 
sional case. In the case of infinitely many degrees of freedom it was 
introduced by Segal [H] (see also [15]). One of the most famous appli- 
cations of this representation we find in the book [3] by Berezin, where 
it has the interesting counterpart to the Fermion case, i.e., canonical 
anti-commutation relations. The canonical anti-commutation relation 
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has a representation which is similar to the holomorphic representation 
where however the field of complex numbers is replaced by the Grass- 
mann algebra. Such a representation is now popular under the name 
of Berezin calculus. Another important application of the holomorhic 
representation is given in the book of Faddeev and Slanov [8] . 

4.1. The holomorphic (Bargmann) representation of the CCR. 

We develop the representation of canonical commutation relation (CCR) 
called holomorphic representation, which seems to be quite useful for 
quantum optics and quantum teleportation of continuous variables. 



The operator of multiplication q and the differentiation p = —id/dq in 
L 2 (IR) = L 2 (R, dq) satisfies the commutation relation: 

[q,p] = qp — pq = il 

on a suitable subspace of L 2 (M). Introduce the operators 

a = (q + d/dq) /V2 = (q + ip) /V2, a ] = (q - ip) 

then 

[a, a f ] = 1, q = (a + a f ) /V2, p = (a - a ] ) /y/2i. 
This representation of the commutation relation is called the Schrddinger 
representation. The function f(q) = e~ q I 2 is a solution of the equation 

= V2af(q) = (q + d/dq)f(q). 

We consider the space L 2 (IR, e~ q dq/ ^n) and the unitary operator U 

U : L 2 (R,e- q2 dq/V2^) 3 g(q) -> ^Y^gisi)^ 12 G L 2 (R,dq). 
Then we have 

iPqU = q, U\d/dq)U = d/dq - q, 

b = U ] aU = U\q + d/dq)U/V2 = 2- 1/2 d/dq, 

&t = u ] a ] U = U\q - d/dq)U/V2 = 2~ 1/2 (2g - d/dq). 

This representation L 2 (R, e^g/v 7 ^), b = 2~ 1 / 2 d/dq, = 2~ 1 / 2 (2g- 
d/dq) is called the modified Schrddinger representation. 
Next we consider 

u = x + iy, L 2 (C,^), ^ = e-^ = e-^ 2 )^, f dp, = 1, 

27T1 7T J £ 

and its subspace Ti generated by holomorphic functions of u (anti- 
holomorphic functions). For /, g 6 TC define the inner product by 

dudu 



(f\g) = j f(u)g(u)e ^ 

Then the multiplication operator f(u) uf(u) and the differential 
operator f{u) — > d/duf(u) are adjoint to each other: 

dudu f „ , „ r „-,., —, r ^ , x dudu 



f(u){d/dug(u)}e- uu — = - / d/du{e- uu f(u)}g(u)- . 

Z7T1 / Z7T1 



/ 
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,dudu f — -r-^ , , _^,dudu 



f{u)ug{u)e- au< ^ = J uf(u)g(u)e 



27Ti ' 



where we used the relation 8/ duf{u) = because f{u) is a holomorphic 
function of u. If we put 

a = d/du = (l/2)(d/dx + id/dy), = u, 

then a and a) satisfy the commutation relation [a, at] = 1. We call this 
representation the holomorphic representation. In this representation, 
{u n /Vn\}™ =0 is an orthonormal basis. In fact, let m < n. Then we 
have 

u m u n e _ Uu dudu r u rn n-l) n d n \ Uu dudu 



\fm\ VnJ. 27ri J y^T y \fn\du a J 27ri 



d n u m \ _ r „„dudu 

g aw — o, 



VnJ.du n VmJ.J 2ix\ 



m il ■ 



and if m > n, 



u m u n _ Uu dudu I* / d m u n \ _ Qu dudu 



'm\\fn\ 27ri J \^/^\du m \fn\ ) 2vri 

Therefore, if /(iZ) = X^=o a ™^™ tnen H/ll 2 = X^=o Kl 2n! - Tne co- 
herent state is an eigen state of the annihilation operator. Also in the 
holomorphic representation, 

^f(u) = «/(«), /(«) = Ce^ = CJ2 

n=0 

oo I |2 n 

C = e-H 2 /2 ^ HJII2 = | C |2 £ N_ = | C |2 e |aP = L 

n=0 n ' 

Hence, / is normalized by choosing C = e~' a ' 2 / 2 . Furthermore, in the 
holomorphic representation, the integral kernel K(u, v) of the identity 
operator is 

K(u,v) = e av . 

In fact, for a G C, 

e *v-n e -^dvdv = f ^ (a^^dvdv = qU 



2m / ^ k\ 2vri 

k=0 



implies 
(4.1) 



/ e™f(v)e-^ = I e-^a^e-|^ = £a n c*" = f(a). 

Because of symmetry in v and v, one also has 
(4.2) f e»e"*^ = f(a). 
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The equality 

a) m a n f(u) = u m -^- [ e Uv f(v)e-* vd ^ = f u m v n e Uv f(v)e^ vdvdV 



du n J 2ni J 2n\ 

shows that the integral kernel of the normal ordered monomial o) m a n 
of a* and a is u m v n e uv . 

The holomorphic representation for the case of n degrees of freedom is 
briefly described below: 

The variables, the Hilbert space and the measure now are: 



u={u u ...,u n ), L 2 (C\o> n ), oX = n e ~^^T' 

3=1 



4.2. Integral kernels of basic operations in the holomorhic rep- 
resentation. For the mathematical description of the devices for the 
manipulation of photon states we need several results for the kernels 
of various linear transformations of the basic creation and annihila- 
tion operators. The relevant results are stated in this subsection. The 
proofs are contained in the appendix. 

These results are special cases of results given in [3] where the infinite 
dimensional counter part is given. The proofs of these results in [3] are 
based on the use of functional integration. Though the basic strategy is 
the same as in [3], the proofs given in our appendix only uses elementary 
mathematical tools so that these result become more easily accessible. 

Theorem 4.1. The linear canonical transformation 
(4.3) b, = a, + ./•, b) = a) + 

is implemented by the unitary operator U whose integral kernel U(u, v) 
is 

n ( n ~\ 

U(u,v) = cexp^(%u i +T; i / i -%/ i ), c = 6exp I (-1/2) fofo > , \6\ = 1. 

3=1 I 3=1 ) 

Theorem 4.2. The linear canonical transformation 

n n 
k=l k=l 

is implemented by the unitary operator U whose integral kernel U(u, v) 
is given by 

U(u,v) = cexp Q(« u) ( ^ ^2 ) ( I )) , A h = 'A k K 
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Theorem 4.3. Let fj be complex numbers and 

n 

i=i 

be a self-adjoint operator. Then the kernel function U(u,v) of e ltH is 
given by 
(4.4) 



U(u, v) = cexp 



n f 1 " - 1 

3=1 I j=l ) 



Theorem 4.4. Let B and C be n x n matrices, and 

H = i(aW + aBa + 2a ] Ca) 
be a self-adjoint operator. Further, denote 

A 









f $ 


¥ \ 




-ft 


, e^ = [ 


v * 





Then the kernel function U(u, v) of e ltH is given by 

(4.5) U(u,v) = cexp Q(„ u) ( f 2X a» ) ( S )) ' ^ = ^ 
with 

(4.6) A 22 = A 21 = S" 1 , A 11 = M" 1 , c = (det $e itc, )- 1/2 . 
5. Photon states and devices for manipulating photon 

STATES 

We consider a self-adjoint operator 

Hu ser = i(aa - aa)) 

and an operator D(a) = e lHlasei ^ a ' called the displacement operator 
which generates from the vacuum state a state called the coherent state 
which is considered to represent the laser-beam state. In order to have 
the kernel U(u,v) of the operator e lHlascr ^ 9 \ we use Theorem 14.31 i.e., 
n = 1 and / = — ia. The kernel U(u,v) of the operator D(oe) therefore 
is 

U (u, v) = e~'°' / 2 exp{uv — ctv + cm}. 
In the holomorphic representation, its action on a state can be calcu- 
lated explicitly 

(D(a)f)(u)= j U(u,v)f(v)e-*^ 

dvdv 



e-^jeMuv-a V + au}me 
:e -M72 e «i f ew {(u-a)v}f{v)e- 



2vri 

,dvdv 
2vri 
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(5.1) =e~ H2/2 e aii f(u-a). 

If f(u) = 1 which corresponds to the vacuum |0), then 

\ a ) = D(a)\0) = e - H/2 e aa 
is the coherent state and satisfies the following completeness relation: 

(5.2) f \a)(a\^ = ( e^W^ = e*, 



2vri J 2vri 

where we used the fact that e uv is the kernel of the identity operator I. 
Next consider a self-adjoint operator H parSu (g) = ig(a 2 — a^ 2 ), the gen- 
erator of parametric amplification. In order to have the kernel U(u, v) 
of the operator e lHpara ^, we use Theorem 14.41 i.e., B = —ig and C — 0, 
and 



o j ' ^ * $ y F y \^ - i o 

,/l0\ . , x f -1 \ / cosh to -sinht# 
= cosht^ 1 j+smht^ _ 1 Q J = ^_ sinh ^ cosh ^ 

A 22 = = t&nhtg, A 21 = = cosh -1 ^, A 11 = = -tsahtg. 

The kernel U(u, v) of the operator e lHpara ^ therefore is 

U(u,v) = (coshg) _1//2 exp(l/2){tanh gu 2 + 2 cosh -1 #w — tanhgt> 2 }. 

In the holomorphic representation, its action on the vacuum state can 
be calculated explicitly 

= (cosh #)~ 1/2 exp( 1/2) {tanh gw 2 } 



x J exp(l/2){2 cosh 1 guv — tanh gv 2 }e 



2vri 



(5.3) = (1 - tanh 2 g) 1/4 exp( 1/2) {tanh gu 2 } 

where we used the relation (14. 2p . The state (15. 3p is called the squeezed 
vacuum with squeezing parameter g. 

Let us consider the beam splitter. The generator of the beam splitter 
H\> s (9) is defined by 

Hbs (0) = i9(a\a 2 - aia\). 

Then applying the Theorem 14.41 again, we get the kernel U(u, v) of 
e lHhs w as follows 
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® JtA _ ( C~ iC 

$ J \ e itc 

e -itc = e wc = expt0 ^ ^ J^=cost0^J J^+sinte^ _ x Q 

/ costO smt9\ 
\ -sintO costO J ' 



1 22 



A 11 = 0, A 21 C ° S ? " Sin / V c = (det Se^-V* = i. 

y sinf cosft / 



The beam splitter for 6? = 7r/4 is called the half-beam splitter whose 
generator is i^hbs = i(7r/4) (a{a 2 — a\a 2 ). The half-beam splitter trans- 
forms the state f(ui,u 2 ) into 

U(u,v)f(v u v 2 )e^f[^ 



I 



y"exp-^={(Mi + u 2 )v 1 + (-Mi + u 2 )v 2 }f(vi,v 2 )e- vv Y[ Vj "° 3 



2ni 

3=1 



(5.4) =/((m 1 +m 2 )/v / 2,(-mi + m 2 )/v / 2). 



5.1. Balanced homodyne detection. We calculate the mean value 
of the difference of photon beams N\ — N 2 for the state 

e -iH hbs l^^l^^ the state after passing through the half-beam splitter. 
Since 

e iHhbs a ie - iHhbs = (ai + a 2 )/V2, e iHhbs a 2 e~ iHhbs = (-ai + a 2 )/V2, 
(H ® (^le^^i - A^ 2 )e- i ^|V'i> <8> ^2) 
= (ipi\ ® (^2\(a\a 2 + ai4)|^i) ® |^)- 
If \ip 2 ) is a coherent state \a 2 ), then this identity is continued by 

= <-0i| ® (a 2 |(c4a 2 + ai4)|^i) <8> |a 2 ) = (^1 1 all ^1)^2 + (^iKI^i)^ 
= (M<7i - ipJ/yfiMlatle** + (Vil^i + ipiVv^OHe-"* 2 
= a/2 <-0i I Q'x I "0x> I QJ2 1 cos# 2 + V2(ipi\pi\il)i}\a 2 \ sin6 2 
= I (Alltel cos6» 2 +p! sin 6*2)1^1). 

Thus the balanced homodyne detection measures v / 2|« 2 | times of q\ 
and p 2 according to the phase of the coherent beam a 2 . 
The displacement operator D(— a) is realized by modulators and the 
beam splitter. For \<f>) e H 2 and \j3) G H3, we send |0) ® \f3) through 
a beam splitter B(6) = e iH "^ e) . Let 

|0)®| / 3) = 0(M 2 ) e -l' 3 l 2 /V ii3 . 
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Then 

B{6){\<f>) <g> = B{6)<P{u 2 )e-^ 2l2 e^ 
= (j)(u2 cos 9 + u 3 sin 0) e - |/3|2/2 e /3( -" 2 sine+ " 3 cose) 
= 0(m 2 cos 9 + u 3 sin #) e - |/3|2/2 e-^ 2 sine e ^ 3 cosS . 

(I®(/5|)fl(W>®l/?> 
V W3 0(«2 cos + u 3 sin ^e"^ 2 sin9 e^ cos V* 3 " 3 

= e -|/9|» e -^ aB tofl y e /3n 3e ^3cos^ ( - 2cos ^ + - 3sin ^*^3 

= e-We-P** sine e^ cose <P(u 2 cos 9 + sin 9) . 
Let # — > and /3 — > oo such that (3 sin 9 — > a. Then 

|/3| 2 (1 - cos6>) -> |a| 2 /2 

and 

e -|/3| 2 e -/3« 2 smfl^cose^^ cos + ^ s i n 0) 

-> e- |a|/ V aS2 0(M 2 + a) = D(a)\4>). 

6. Teleportation of continuous quantum variables 

The proposal for continuous variable quantum teleportation was first 
made by Vaidman [T7] in 1994, and then by Braunstein and Kimble 
[6] in 1998, and experimentally demonstrated by the Caltech group, 
Furusawa et al. [10] in 1998. Now we present the teleportation of 
continuous quantum variables in a parallel way as the teleportation of 
qubits in Section 2. The numbering (1), . . . , (4) corresponds to that in 
Section 2. Figure 2 is similar to the Figure 1 and these figures illustrate 
the correspondence of the two cases of teleportation. 

(1) First, an entangled state as a counterpart to an EPR pair is 
produced. We prepare such a state using parametric amplifica- 
tion and a beam splitter. By parametric amplification we create 
a pair of squeezed vacuum states U paxa ,(g)\0) = e lHpara ^^|0) and 
f^para ( _ flO|0). In the holomorphic representation one has (see 

(l-g 2 )- 1/4 f/pa r a(^)|0) = e^/ 2 , 

(1 - <7T 1/4 £W-2)|0) = e"^ /2 , q = tanhs. 

Then we send these states through the half-beam splitter (see 
f)5.4p ). Now Alice and Bob share the state {u\ is Alice's variable 
and U2 is Bob's variable) 

|tf ) = e 9(«l+«2) 2 /4 e -< ? (-«i+« 2 ) 2 /4 = e ,»ii 2) q = tanhg _ 
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Victor 



Bell state measurement 



7^ 







e- Q l— -y 





p+ 



Alice 



Bob 



Quantum entanglement 



D(x_ + ip + ) -- V»(U2) 



Unitary transformation 



Measurement 



D(a) 



Displacement operator 



Half-beamsplitter 



Figure 2. Quantum circuit of continuous variable teleportation 



Victor gives Alice a state ip(uo) to send to Bob. The state of 
the total system is 

oo 

= ^{u )e qai ^ = |^)o ® ^Q n \n)i ® \n) 2 , q = tanh^. 

n=0 

This state corresponds to the state 

|^o> = ^>v®^(|0>a|0) b + |1)a|1)b) 

of the qubit case. The generalized state (for some background 
information on generalized states see the first part of the Ap- 
pendix) 



= -- 1/2 E i! ^r = ^ 1/2 EHo^i^ 



n=0 



n=0 



corresponds to the Bell state \/3 Q0 ) = 2 x l 2 J2n=o \ n )o ® \ n )i of 
Section 2. The other Bell states \f3ij) correspond to the states 



J2(D(a)\n) ) ® |n)i = e -H 2 /^ e ^0-a)v le av = e -\a\y2 e v v le -av le 



n=0 
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(see (EH))- 

In the Appendix we show that 

oo 

(6.1) {tt- 1/2 ^T(D(a)\n) ) <g> |n)i; a G C}, a = x_ + ip+ 

n=0 

is the generalized Bell basis, i.e., a complete orthonormal system 
in our Hilbert space. 

Then we send this state through the half-beam splitter real- 
ized as the unitary operator e tHhbs , and obtain (see (15.4j) ) 

e -(o : i+^)/2 e -ug/2 e u?/2 e ^_«o e i^P+«l = n V2\ x _) g, 

where 

is the generalized eigen-state of the operator 

1 / d 

with eigen- value X-, and 

\ p+ ) = 7r -l/4 e -p 2 f /2 eS ?/2 e ^p + ^ 1 

is the generalized eigen-state of the operator 

d 



Pl V2i \dui Ul 

with eigen- value p + . In the Appendix the orthogonality rela- 
tions 

(6.2) (x'_\x-) = 5(x'_ - x_), (p' + \p+) = 5(p' + -p+). 
are shown. Thus we have 

oo 

(6.3) e lHhbs ^2(D(a)\n) ) ® \n)i = n 1/2 \x^) ® 

n=0 

and 

oo oo 

o(^IP(«) t ® i(m|)(^(L>( a )|n) ) ® |n)x) 

m=0 n=0 

= O-l ® (pVI): (M ® = 5{x'_-xJ)5{p\-p + ). 
The state | , 0}y ® |Ajo)ab corresponds to 

oo 

|$i) = vr- 1/ V)o ®^q n \n)i <8> |n) 2 . 

n=0 

The relation 

/oo oo 
cfe_dp + ^(D(a)|m)o) ® |m)i (^|£ ) (tt) t ® i{k\ 

m=0 fc=0 
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implies for a = x_ + ip + that 

/oo 



m=0 



;££o<A:|D(a)t|^?M*|n>i|n> a 



fe=0 ra=0 



cfe_c/p + ^(D(a)|m)o) ® ® ^g n (n| J D(a) t |V')o|n}2. 

m=0 n=0 



This is the expansion of with respect to the generalized 
Bell basis (JB1J) . 

(2) Then Alice performs the (generalized) Bell-state measurement, 
and when some Bell state 7r -1 / 2 X]^=o(-^'( Q! )I 77 ')o) ® l n H * s cno " 
sen, then the total system is reduced to 

oo oo 

= TT" 1 / 2 ^(£>(a)|m> ) ® |m>! ® ^g n o(n|D(a)t|V>o|n> 2 . 

m=0 n=0 

If g = 1, the above state is 

oo 

J2(D(a)\m) ) ® |m)i ® D(a) f |^) 2 . 

m=0 

Thus the Bob's state is reduced to £)(a)'| , 0)2 though he can 
not know this. 

(3) Alice sends the classical information a = X- +ip+ to Bob. Then 
Bob sends Bob's state D(aYip through D(a) according to the 
result a = x_ + ip + , obtaining ^ out = D{a)D{a)^ij) = ip- Alice 
sends two real numbers (x_,p + ) to Bob, and Bob gets a state 
ip = Y2™=o a n\ n ) m an infinite dimensional Hilbert space. 

(4) The Bell-state measurement (2) is done by sending the Bell 
basis through the half-beam splitter obtaining the canonical 
basis {\x-) ® \p + );a = X- + ip + G C} of (16. 3D and by the 
measurement of xq and p\ using balanced homodyne detection. 

Since we cannot generate an EPR pair e" 1 " 2 with an infi- 
nite squeezing parameter g = oo, the ideal q — 1 case of the 
squeezed state e quiU2 , q = tanhg < 1, we cannot have com- 
plete teleportation ip out = ip(— ip- m ). We have to measure the 
distance between ip and ipout- For density operators p (positive 
operator with tr p = 1) and a, the fidelity F(p, a), given by 



F(p, a) = tr v / P 1/2 ap 1 /2 
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is an accepted measure for such a distance. Let \if)) = \j) be a 
coherent state, and 



J2q n D(a)\n)(n\D(ay\j). 



n=0 



For p = | {tI an d cr 



VpV^pV! = = Vl7>(7l0>(0l7>(7l = I<7l0>li7><7l 

and the fidelity F(p,a) equals |(7|0)|. The inner product (7|</>) 

is easily calculated, using relation (14. 1J) : 

oo i> i> 

,|2/n _|„,|2 



d\ ( i>) = J2i n ^\ D ( a )\ n )( n \ D ( a )^)= / / 



e ~\ a \ l 2 e ~ au e -\l\ /2 e ~f(u+a) e quv 



xc - M y2 c - a v c -h\V2 cl (v + a) c -uu c -vy dududvdv 

2vri 2vri 

= exp(l — g){«7 + «7 — |a| 2 — |7| 2 }. 

7. Experiment, Controversy and Locality 

7.1. The Experiment. Here, we briefly explain the experiment by 
the Caltech group, Furusawa et al. [ID] . The teleported state is not an 
optical beam itself but a modulation sideband of a bright optical beam 
generated by an electro-optical modulator, because the frequency of 
the optical beam is too high (uj/tt = 300THz, wavelength 1000 nm) to 
handle directly. So, the optical beam is treated only as a carrier and 
the quantum states are discussed using sideband frequency. The light 
from a single-frequency titanium sapphire (TiAl 2 03) laser at 860 nm 
(frequency u£) serves as the primary source for all fields in the exper- 
iment. 90 percent of the laser output with frequency ul is directed to 
a frequency-doubling cavity to generate blue light at 2u>l- This output 
then splits into two beams that serve as harmonic pumps for para- 
metric down- conversion, 2ujl — > ojl ± Q, within the optical parametric 
oscillator (OPO). 

(1) Thus the two independent squeezed beams which are supposed 

-2 -2 

to be represented by e qUl and e _<? " 2 at 860 nm are generated 
by optical parametric oscillators (amplifiers), where q = tanhg 
and g is called the squeezing parameter. These beams are sent 
through a beamsplitter obtaining the EPR beam e quiU2 and the 
outcomes are sent to Alice (variable u\) and Bob (variable 7Z2). 

Victor generates a coherent sideband \tp) = e -\ a \/ 2 e au » at 
frequency (00 ± Q)/2rr with Q/2tt = 2.9 MHz by an electro- 
optical modulator, and sends it to Alice. 

If e ltu) represents the carrier and ae^ the modulation, then 
the total system is represented by 

ae ±itn e ito = ae i{ul±n) . 
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M x | M p | optical modulator QPQ optical parametric oscillator 

LO (local oscillator) 



input 




output 



Balanced 
homodyne detector 



-p Half-beam splitter 



99/1 beam splitter 



FIGURE 3. Experimental scheme for quantum telepor- 
tation performed by Caltech group. The laser field is 
shared by all parties. 



(2) (a) Alice sends two beams, e~' a ' //2 e Q "° and her EPR beam, 
through a half-beamsplitter e lHhhs discussed in Section 5, 
obtaining 



e ~\a\/2 e o(M +ui)v / 2 e <j(-Mo+Mi)«2 / V2 



20 



S. NAGAMACHI AND E. BRUNING 



(b) Alice uses two sets of balanced homodyne detectors (D xo , D pi ) 
discussed in Section 5 to make a joint measurement of the 
amplitude x and pi discussed in Section 6. 

(3) Alice sends the outcomes (x_,p + ) of the detectors to Bob. 

(4) After receiving this classical information from Alice, Bob is able 
to construct the teleported state p out . Bob generates a sideband 
beam at frequency (a; ±Q)/2n by two electro-optical modula- 
tors M x (amplitude modulator) and M p (phase modulator) with 
suitable complex amplitude \(3) = e'^^e^ 3 . Then he sends |/3) 
and his EPR beam through a beamsplitter of refractivity 0.99 
obtaining the state p out . 

(5) Victor detects the state p out by his own balanced homodyne 
detector Dy and compares it with his original state 

The result obtained by Furusawa, A., et al. is that the fidelity F(\ip), p ont ) 
of the states \ip) and p ont is 

F(\ip), Pout) =0.58 ±0.02. 

Bowen, W.P., et al. [5] had the fidelity of 0.64 ± 0.02 and Takei, N., et 
al. p] achieved the fidelity of 0.70 ± 0.02. 

7.2. Controversy. After the experimental demonstration of continu- 
ous variable quantum teleportation (CVQT) was reported [TP] , there 
was a controversy over its validity on the ground of intrinsic phase inde- 
terminacy of the laser field [T2] . The laser field is often assumed to be 
a coherent state having a fixed phase, but [12] shows that the steady- 
state solution of the master equation in the quantum theory of the laser 
shows that the phase of the laser field inside the cavity is completely 
unknown and genuine CVQT cannot be achieved using conventional 
lase sources, due to an absence of optical coherence. Furthermore, the 
same laser source is used for (i) producing Victor's state for telepor- 
tation, (ii) pumping the nonlinear crystal that produces a two-mode 
squeezed light field serving as the shared EPR state, (iii) supplying lo- 
cal oscillator (LO) fields for both of Alice's homodyne measurements, 
and (iv) providing Bob with a coherent field to mix with his portion of 
the EPR beam to reconstruct Victor's state. 

The ideal scheme for CVQT [6] is explained in Fig. 2 and an explana- 
tion of the experimental simulation of CVQT [10] is provided in Fig. 
3. 

For the claim [12], there appeared a counterargument [IB] which says 
that the standard description of the laser field used in [12] is insuffi- 
cient to understand CVQT with a laser. They found that the laser 
light has the random phase only in the cavity, and outside the cavity, 
it has similar phase character as the coherent state has, and therefore 
a conventional laser can be used for CVQT. 

Furthermore, [9] says that the laser field outside the cavity is a mixed 
state whose phase is completely unknown, but CVQT with a laser is 
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valid only if the unknown phase of the laser field is shared among 
sender's LOs, the EPR state, and receiver's LO. 

7.3. Locality. Let fi(x) be two functions square integrable functions 
on Euclidean space, i.e., elements in L 2 (M 3 ), and suppose that the 
support of the function fi is contained in a bounded set 0« of M 3 . 
Then for a system in the state represented by the function fi(x), the 
observable x (the position of a particle) is always observed in the set 
Oj. In such situation, we often say that the state fi is localized in Oj. 
The state of two particles at x\ and £2 which are separated by a long 
distance is expressed by the tensor product f{x\,x-i) = fi{x\)f2{x2) of 
two functions fi (i = 1, 2) whose supports Oj (i = 1, 2) are separated 
by a long distance. 

Ideal quantum teleportation can be considered as the following pro- 
cess: The information contained in a quantum state |^)i n (which is 
localized in a region 0{) is sent to another region O2 separated from 
0\ by a long distance and produces a state iV^out (which is localized 
in O2) and contains all the information of iV^in ■ 

In the framework of teleportation of qubits in Section 2 where a finite 
dimensional Hilbert space is used, it seems difficult to define a state lo- 
calized in some bounded region. In the framework of CVQT in Section 
6, though an infinite dimensional Hilbert space is used, we manipu- 
late only monochromatic laser beams. When the frequency of the laser 
beam is determined, the location of the laser beam is completely un- 
known. Up to now, one has considered quantum teleportation without 
the notion of locality. In a precise formulation of quantum teleporta- 
tion which takes locality into account, a quantum theory of infinitely 
many degrees of freedom might be needed, e.g., a quantum field theory 
might be necessary. 

It seems that the theory of lasers is not satisfactory in its appli- 
cation to quantum teleportation, and therefore, many controversies 
appear around CVQT which uses laser beams. In this context to, a 
theory of quantum teleportation based on quantum field theory might 
be preferable. 

Appendix A. 

A.l. Generalized vectors and generalized eigen-states. The the- 
ory of the generalized vectors is developed in [TH]. Let 5(IR) the 
Schwartz space of fast decaying infinitely often differentiable functions 
on the real line R and L 2 (R) the Hilbert space of square integrable 
functions on R. Then 

<S(R) C L 2 (R) C «S'(R) 

is a rigged Hilbert space (see [T9j ) . Here iS'(R) denotes the topological 
dualofS(R). 
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The elements F of S' are called generalized vectors. Suppose that an 
operator A on L 2 (R) maps <S(R) into «S(R). Then a generalized vector 
F x G <S'(R) is called a generalized eigen- vector of A corresponding to 
an eigen- value A, if 

F\{A<p) = AF A (0) 

holds for every G <S(R). 
Let 

E x = {FeS'{R);F{A^ = \F{ ( p)} 
be the space of generalized eigen-vectors of A corresponding to the 
eigen- value A. We associate with each element G <S(R) and each 
number A a linear functional A on E\, taking the value F\((f>) on the 
element F\ of E\. We call the correspondence — > A the spectral 
decomposition of the element corresponding to the operator A. 
If A = implies = 0, then we say that the set of generalized eigen- 
vectors of the operator A is complete. 

Let p = —id/dx be the self-adjoint generator of translations in L 2 (R) 
(more accurately p is the self-adjoint realization of the differential op- 
erator —id/dx on <S(R)). Then the identity 

{-id/dx)e lXx = Xe iXx , A G R 

shows that e lXx is an eigen-function of the operator p corresponding 
to the eigen- value A. The function e lXx does not belong to L 2 (R), but 
belongs to S'(R) 3 F x = e iXx . The relation 



F\(p(/>) = / e lXx (— id/dx) (J)(x)dx = / {—id / dx)e lXx <\>{x)dx 



A J e iXx (p{x)dx = AF A (0) 



shows that F\ = e tXx is a generalized eigen- vector of p corresponding 
to the eigen-value A. The spectral decomposition A of is the Fourier 
transformation of 0. 



A = F A (0) = J e iXx <p{x)dx = I e~ iXx <f){x)dx = 0(A). 
The Fourier inversion formula 



m = ^ / e iXx 4>{\)d\ 



shows that A = implies = 0, that is, {e lXx ; A G R} is a complete 
set of generalized eigen-vectors of p. 

Though it is not possible to find the value of the momentum p (con- 
tinuous variable) to be precisely A, it is convenient to say that if the 
value A of the momentum p is found for the state |0) = 0(x), the state 
after the measurement is |A)(A|0) = e lXx (p{\) /2-k in the same way as 
the discrete variables, where |A) = e lXx /V2n is a generalized vector. 
Actually, we can only say that A is contained in the interval [a, b]. In 
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that case, using the above generalized eigen-states, the state after the 
measurement is 



^A|A)(A|0)/y^rfA|(A|0)| 2 

. Thus we have the so-called projective measurement, whose name 
comes from the fact that P = J dX\X)(X\ is a projection. If we think 

that the edge of [a, b] is not sharp, then we take M = dXx(X) | A) (A| as 
a measurement operator, where x{X) is a C°°-function with the support 
contained in [a — e, b + e] and x(A) = 1 for A G [a + e, b — e] and e > 0. 
Then we have the following generalized measurement (see [20J). 
The state after the measurement is 



M|0)/v/(0|M*M|0>. 



If we introduce the operators 



y/2\ dx J ' a " y/2 V dx) ' 



the topology of S (R) is defined by the system of norms 

||0|| r 2 = (0|(l + a*ar|0> 
for any r G N (see Theorem V.13 of |21j). 

In our case of a rigged Hilbert space $ C 7i C $ is the set of vectors 
of the form 



n=0 



satisfying 



n=0 

for any r G N. $ is a countably normed space. Since 

oo 

^(l + n 2 ) r |A|7n! < oo, 



n=0 



e G $ for any A G C, the complex numbers. 
For \q\ < 1 we have 



oo 

,2n 



j>»o®Hiii 2 =;>> 2 W 2 

n=0 n=0 ^ 



but for q = 1 this series is divergent: 

oo oo 

ii ^2 ® i n ^n 2 = 1 = °°- 



n=0 n=0 
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In this sense the vector 



-l/ 2 ^o«i 



7r ->-e = ^rm =ir i/»r Ho^ln)! 
^— ' n! 

n=0 n=0 

is a generalized vector. The unitary operator e lHhhs sends the general- 
ized vector 

oo 

Y^(D(-a)\n) ) <g> |n)i = eH^/V^e^e^ 

ra=0 

to the generalized vector 

= e -(* 2 _M)/ 2 e -^/ 2 e *?/ 2 e ^-«o e ^+«i = TT 1 / 2 ^-) <g> 

where a = x_ + ip + and 

= ^-i/4 e -*l/2 e -«§/2 e >/2*_«o 

is the generalized eigen-state of the operator 

1 / d 

with eigen- value x_, and 

= 7r -l/4 e -p 2 + /2 eS ?/2 e iV2 P+ « 1 

is the generalized eigen-state of the operator 

1 / d 



with eigen- value p+. In the usual Schrodinger representation, \p + ) is 
represented by {2n)~ l l 2 e ip+x and is a generalized eigen-vector of p — 
—id/dx as explained above. Since Ti is a separable Hilbert space, has 
a countable basis. But for the completeness, 7i must have uncountably 
many orthogonal generalized vectors. 

A. 2. Proofs of theorems 1-4, Section 4. 

A. 2.1. Proof of Theorem 1. First, consider the transformation 

h 3 = "l + //• h ) = a ) + fj 

for complex numbers fj. We want to have a unitary operator U such 
that 

bjU = Uaj, b)U = Ua) bj = U aj U\ b) = Ua]rf. 

We assume that the operator U is defined by a kernel U(u,v), that is, 



„ n 

(Ug)(u) = J U(u,v)g(v)l[e-°™ 



dvjdvj 



2m 

3=1 
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Then we have, by integration by part, 



2m 



3 \ j=l ) ' j=l 

and similarly 



. dvjdvj 
2ni 

r- 



/IV 
U(u,v)v jg (v)l[e-^ 
j'=i 



Thus, it is clear that one should have the following correspondence: 

d 

CLjU <-> — U(u,v), CLjU <-> UjU(u,v), 
Uj 

d 

Ucij <-> VjU(u,v), Ua\ <-> — — U(u,v). 

J OVj 

If we assume that the kernel of the operator [/ has the form 

n 

[/ (w, f ) = cexp ^^(fZji'j + f + Mj^j) 

3=1 

for complex numbers 0^ and ipj, then we get indeed 



bjU (u, v) = ( + fj ) U (u, v) = VjU(u, v), 



d 

b*U{u,v) = (Uj + fj)U(u,v) = —U(u,v). 

Since 



d 

U(u,v) = [Vj + ipj)U(u,v), 



duj 



d 

—U (u, v) = (uj + (pj)U(u, v) 

equating the coefficients for Uj and Vj we find ipj = —fj and <f>j = fj. 
Thus the result is 

n 

U (u, v) = cexp J2(ujVj + Vjfj - Ujfj). 

3=1 
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Applying U to the vacuum f(v) — 1 produces 

/n n 
c exp Y^*M + vjfj - Ujfj) J] e- "* 



dvjdvj 



2m 



= cexp<^ - ^Ujfj 



n 

3=1 



r ( n 1 n 

ii 2 = | C | 2 y exp|-x)[ti^+ Vi]|n e " flJU - 

|c| 2 y exp - J ^[(iij + fj)(uj + fj) - fjf. 



n=l 



The normalization of F, i.e., ||F|| = 1, requires 

c = 0exp j(-l/2)£/^-j, \9\ = 1, 
and we have the following theorem. 



A. 2. 2. Proof of Theorem 2. Now we consider more general linear canon- 
ical transformation (Bogoliubov transformation) 

n n 

b J = ^2^jka k + y jk al), b* = ^(Qjkal + * jk a k ) 

k=l k=l 

satisfying 

[bj,b k ] = [bj,b k } = 0, [bj,bl] = 5 jk , 

and find the kernel U(u, v) of the unitary operator U which implements 
the linear canonical transformation 

bjU = Ua,j, b*U = Ua*. 

From the above commutation relations, we have 

n 

= [bj, b k ] = [®H a i + *i' a *> ®k m a m + 

i=l,m=l 

n 

= fe*ta[fli,a;]+¥tak>m]) 

l=l,m=l 
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^ {^jl^km^lm ^ jl^km^lm) ^ i^jm^km ^jm^km) 
/=l,m=l m=l 
n 

5 jk = [bj, b* k ] = ^ [®ii a i + ^H a *n ®k m a* m + ^ km a m ] 



i=l,m=l 

n 



= Mtah.aJ + ^K.aJ) 

i=l,m=l 

^ ^ {^jl^km^lm ^ jl^ km&lrri) ^ ^ (^jm^fcm ^ jm^ km) ■ 
i=l,m=l m=l 

These relations can be described by the matrix notation: 
(A.l) = - / = $$*- tftf*. 

We want to have a unitary operator [/ such that 

bjU = Uaj, b*U = Ua*. 
We assume that the kernel of the operator U has the form 

U(u,v) = cexp (hv u)( A ^ A* = M*. 



Then we have 

n 

du k 



d 

bjU(u, v) = ^2(®jk-^— + ^ jk u k )U(u, v) = vjU (u, v) 
k=i 



Since 

d 



du k 



U(u,v) = J2( A l m v m + Al 2 m u)U(u,v), 



m=l 

equating the coefficients for ui and vi we have 

$A 22 + ^ = 0, $A 21 = /. 
In the same way, from the equation 

b*U(u,v) = ^2($ jk a* k + V jk a k )U(u,v) = —U(u,v), 



k=i 



j 



we deduce 

q,A 22 + § = A 12 , ^A 21 = A U . 

This gives 

A 22 = A 21 = $~\ A 11 = 

Applying U to the vacuum |0), we have 



m =Ul0) = c J exp (I(„ „) ( £ % ) ( I )) 7^*^ 



2?ri 
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cexp(l/2){-M$^ 1 ^M} / exp(l/2){2m~ 1 v+v^^>~ 1 v-vv} JJ 



3=1 



2m 



The integral is 

exp(-l/2) { (v v) 



/ 



+ 2m~ l v \ J [ 



dvjdvj 
2m 



det 



/ 



-1/2 



x exp 



1 
I 







(l/2)(0 2w$" 1 

where we used the following formula of Gaussian integral: 
-1 



exp 



det 



■(f v) 



A 11 A 12 
A 21 A 22 



+ (A /a) 



J21 ^22 \ 1 - 1 / 2 



A 11 A 12 
Now we calculate 



« 1 )(l/2)^(/ J ./,)( \ 2X \ 22 



3=1 

111 All x _1 



ndvjdvj 
2txi 

fx 
h 



\\F\\ 2 = \c\ 2 / exp-{(l/2)(u$- 1 ^u + u<S>- 1 yu) +uu}Y[ 



dujdu 



i=i 



2tt« 



c| / exp < — - (w m) 



det 



K 
U 

-1/2 



ndujduj 



3=1 



2m 



det 
= det 











(J- 









det 



J 



det 



/ 



= det(J - (l)- 1 ^- 1 ^) = det(J - S- 1 ^*^*- 1 ) 

= det(7 - - I)®*' 1 ) = det($*$) _1 , 

where we have the relation ( 1A.1I) . The constant c is calculated to be 

c = #(det$$T 1/4 , \6\ = 1, 

and we get the following result. 
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A. 2. 3. Proof of Theorem 3. Proof. Consider the operator 
(A.2) aj (t) = e itH a j e- itH , a*(t) = e ttH a*e~ ttH . 

Differentiating with respect to t, we have 

ld(lj(i) \ttU\ „ u\i JtHru „ -\„-ttH 



i dt 



[H,a J (t)] = [H(t),a J (t)] = ^ tH [H,a J }e 



< U "{-fA< m = ~f. 



3 1 



= e itH f 3 e~« H = /,. 
Integrating this system, we have 

a,j(t) = dj — itfj, a*(t) = a* + itfj. 

Thus ( 1A.2I) is a linear canonical transformation of Theorem 14.11 There- 
fore the kernel U(u, v) has the form (14.4j) . In order to find c precisely, 
we differentiate U(u, v) and e ltH and compare them. 



u,v) 



n , - d \ - 

c Yl ( + fodjT. J U ^' v ^ = C Q2^M + fj v i) + itfjfj)U(u, v) 

7=1 ^ 1 = 1 



J 

As a result we obtain 



\j t = itc( J2 and c = exp J -U 2 J2 fjfj | 

3=1 I 3=1 ) 



A. 2. 4. Proof of Theorem 4- Proof. Consider the operator 



(A.3) a 3 {t) = e itH a je - itH , a*{t) = ^ m a*e~ itH . 

Differentiating with respect to t, we have 

ldagl = [H ^ M = [H{tlaM = e *H [H}aj]e -«H 

n n 

= e ltH -{CjkOk + B jk a* k }e- ltH = - £){C> fc (t) + B jk a* k (t)}, 

k=l k=l 

ld^t) = [H ^ a * {t)] = [H(t),a*(t)] = e itH [H,a*]e- itH 

n n 

= e ltH J2i B J^k + C jk a* k }e- m = £){%,a fc (t) + C jk a* k (t)}. 

k=l k=l 
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Integrating this system, we have 

n n 

= ^(^jkO-k + ^jka* k ), a*(t) = y^(#,-jfcQfc + ^fca fc ), 
fc=i fe=i 

where 

(I l) =exp { Zt ( _ ^ ^ )}' 
Thus, Eq. (IA.3[) is a linear canonical transformation. In order to find 
c precisely, we differentiate U (u, v) and e ltH and compare the results. 

Id fide I d /A 11 A™\(v\\ TT( _ 

1 / 9-9 9 \ 

= c _^ + _ B _ +2sc _j £/(s ,„ ) . 

Put v = u = 0. Then we have 

ldc = ly § A 22 = clyAfBu = c-TtA 22 B = -c-Tr^-^B). 
idt 2 ^ 3 v 2 ^ ]l 3 2 2 v ; 

ij ij 

We shall seek c in the form c = (det M) _1//2 . From the formula det M = 

e Tr log M ^ we find 

4- det M = ( ^-Tr log Af j det M = Tr (4- logM ) dot M 
at \dt I \dt 



I I M- 1 ^ ] del .V. 



Thus we obtain the equation for M: 



I Tr (m- 1 ^ =Tr($- 1 ^fi). 



Recalling that the operators $, \I/ are the solution of the equation 

ld_ f $ * \ _ / $ * \ / -C -B 
idt\* ® )~\* ® )\ B C 

we have 

z at 

We now set M = $>e ttc , then we have 

x dM\ / itC r ! fid® itC ^ n jtc 



-Tr M — — = Tr e - lto $ _1 — -e 4to + $Ce 
z \ at / \ \i at 

= Tr (V^" 1 - — e^ + eT^Ce** 7 
\ z dt 

= Tr (e- i4C (-C + $- 1 ^fi)e itc + e^Ce** 3 ) = Tr 

Thus c = (det $e ,iC, )~ 1//2 . Therefore the kernel U(u,v) has the form 

(USD, flU). 
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A. 3. Proofs for Section 6. 

A. 3.1. Generalized Bell basis. First we show completeness: 

n -l f e -^ e ^u le -au le auo e vov le -av le av odx _ dp+ 



;t 1 / e- aa e a ° ai e aivo - ai) e aia °- Vl) e VOVl dx_dp A 



da da 
2vri 



gUOMlg(«0-t'l)(l'0-Ml) e '«Ofl _ g^OfOgWll'l 



This is the kernel of the identity operator (see (14.11) ). Thus (16. ID is a 
complete system. 

In order to show that (16. ip is an orthogonal system, we rewrite the 
vector Y^=o(^( a )\ n )o) ® °f 16.1ft as 

e -| a | 2 / 2 e s oSi e -«Sig«So _ e -\a\ 2 /2 e v V! e -(x--ip + )v! e (x-+ip+)v 

= e ~\ a \ 2 1 2 g^-C^O-Si) e ip+(i)o+^i) _ 
_ g-|a| 2 /2 e -(i'o-Si) 2 /4 e (i'o+?)i) 2 /4 e x_(?}o-i;i) e ip + (ijo+Si) 



A. 3. 2. Proof of Orthogonality relation ( ffi.ffj) . This is a straightforward 
calculation: 

-1/2 -x'f/2 -x 2 _/2 / -^/2_>/2xLuo--fi, 2 ,/2 -^-50.-50^0 ^0^0 



/ x / I^V = ^-1/2^-/2^x1/2 / e -«5/2 e V2xluo e -«S/2 e V2x_«o e 



2m 



n -l/2 e -x'3/2 e -x 2 _/2 f e -2q 2 e V2q(x'_+x-) e iV2p(x'_-x-)^ldP 



7T 



= -1/2 -x'f/2 -xi/2 j -{e-q{x'_+x-))Ap{x'_-x-)<^P 

J 2tt 

= Tt-We^'^-l^eV-^-T'HixL - x-) = 8(x'_ - 
and in a similar way 
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